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ABSTRACT 

Supersonic random motions are observed in dark clouds and are traditionally inter- 
preted as Alfven waves, but the possibility that these motions are super-Alfvenic has 
not been ruled out. 

In this work we report the results of numerical experiments in two opposite regimes; 
M.A ~ 1 and M.a ^> 1, where M.a is the initial Alfvenic Mach number -the ratio of 
the rms velocity to the Alfven speed. Our results show that models with A4a 3> 1 
are consistent with the observed properties of molecular clouds that we have tested - 
statistics of extinction measurements, Zeeman splitting measurements of magnetic field 
strength, line width versus integrated antenna temperature of molecular emission line 
spectra, statistical B-n relation, and scatter in that relation- while models with M.a ~ 1 
have properties that are in conflict with the observations. 

We find that both the density and the magnetic field in molecular clouds may be very 
intermittent. The statistical distributions of magnetic field and gas density are related 
by a power law, with an index that decreases with time in experiments with decaying 
turbulence. After about one dynamical time it stabilizes at B oc n 0A . Magnetically 
dominated cores form early in the evolution, while later on the intermittency in the 
density field wins out, and also cores with weak field can be generated, by mass accretion 
along magnetic field lines. 

Subject headings: turbulence - ISM: kinematics and dynamics- magnetic field 



1. Introduction 

The observation of supersonic motions in molecu- 
lar clouds (eg Zuckerman & Palmer 1974) raised the 
question of how these motions could be supported 
(Norman & Silk 1980; Fleck 1981; Scab & Pumphrey 
1982). Supersonic motions are expected to quickly 
dissipate their energy in highly radiative shocks, be- 
cause of the very short cooling time of molecular 
gas or metal rich atomic gas (Mestel & Spitzer 1956; 
Spitzer 1968; Goldreich & Kwan 1974). 

Strictly related was the issue of the support of 
molecular clouds (MCs) against gravitational col- 
lapse, since it was soon realized that the observed mo- 
tions could not be understood as a gravitational col- 
lapse (Zuckerman & Evans 1974; Morris et al. 1974), 
although MCs contain many Jeans' masses. 

Theoreticians therefore formulated the hypoth- 
esis that MCs were primarily magnetically sup- 
ported (Mestel 1965, Strittmatter 1966; Parker 1973; 
Mouschovias 1976a,b; McKee & Zweibel 1995) and 
interpreted the observed motions as long-wavelength 
hydro-magnetic waves (Arons & Max 1975; Zweibel & 
Josafatsson 1983; Elmegreen 1985, Falgarone & Puget 
1986). It was also shown that the properties of the 
observed 'turbulence' (Larson 1981; Leung, Kutner 
& Mead 1982; Myers 1983; Quiroga 1983; Sanders, 
Scoville & Solomon 1985; Goldsmith & Arquilla 1985; 
Dame ct al. 1986; Falgarone & Perault 1987) could be 
understood if the motions were sub-Alfvenic (Myers 
& Goodman 1988; Mouschovias & Psaltis 1995, Xie 
1997). 

Nevertheless, recent attempts to detect the Zee- 
man effect, in lines of molecules such as OH and CN, 
that probe regions of dense gas (Crutcher et al. 1993; 
Crutcher et al. 1996), resulted in a number of non- 
detections of the effect, and therefore in rather strin- 
gent upper limits for the magnetic field strength, de- 
spite the fact that regions expected to favor detec- 
tions were targeted. It is therefore possible, that the 
field strength in dense molecular gas is weaker than 
assumed in theoretical studies. 

In the present work, we study the dynamics of 
models that share some properties with molecular 
clouds, using numerical solutions of the equations of 
3-D compressible magneto-hydrodynamics (MHD) in 
a regime of highly supersonic random motions. We 
idealize the system by omitting self gravity, ambipolar 
diffusion, radiative energy transfer, etc., and concen- 
trate instead on the aspects that are directly related 



to the random, supersonic nature of the motions. The 
main limitation of our numerical models, compared 
with MCs, is the absence of gravity. We have excluded 
gravity on purpose (our code is capable of handling 
self-gravity) , because one of the aims of our work is to 
show that (magneto-)hydrodynamic processes alone 
are able to explain many of the observed properties 
of MCs. Although gravity is certainly responsible for 
the final collapse of high density regions into stars, 
we suspect that supersonic random motions are re- 
sponsible for setting up many of the properties that 
characterize molecular clouds, with only minor con- 
tributions from gravity (Padoan 1995; Padoan, Jones 
fc Nordlund 1997; P adoan, Nordlund & Jones 1997, 
Padoan ct al., 1998b ). The importance of supersonic 
motions in fragmenting the gas is certainly apparent 
on small scales, where young and probably transient 
clumps are not originated by gravitational instabil- 
ity (Falgarone, Puget & Perault 1992; Langer et al. 
1995). 

Previous numerical studies have shown that com- 
pressible turbulence can qualitatively explain several 
observational properties of MCs, even if the effect of 
gravity and magnetic fields are not included. The first 
two-dimensional (2-D) simulations of turbulence with 
rms Mach number larger than one were performed by 
Passot & Pouquet (1987). These were the first sim- 
ulations where shocks were shown to develop inside 
a turbulent flow. It was immediately recognized that 
shocks might have been responsible for the fragmen- 
tation of the density field inside MCs, and especially 
for the origin of their filamentary structure (Passot, 
Pouquet & Woodward 1988). The importance of 
shock formation inside mildly supersonic flows was 
later confirmed in 3-D simulations by Lee, Lele & 
Moin (1991). Kimura & Tosa (1993) simulated the 
passage of a strong shock through a turbulent molec- 
ular cloud, and found that this process can gener- 
ate dense clumps, with a power law mass spectrum. 
Vazquez-Semadeni (1994) made use of 2-D numerical 
simulations to show that supersonic turbulence gen- 
erates a very intermittent density field, reminiscent 
of the clumpy nature of MCs. The density field was 
also found to be self-similar, which could be the rea- 
son for the hierarchical structure of MCs (Scalo 1985, 
Houlahan & Scalo 1992, Elmegreen & Falgarone 1996, 
Elmegreen 1997). Falgarone et al. (1994), analyz- 
ing the numerical simulation by Porter, Pouquet & 
Woodward (1994), argued that the properties of the 
profiles of molecular emission spectra from MCs can 
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be interpreted as arising from turbulent motions. 

Other numerical works included gravity in the tur- 
bulent flows, yet without describing the magnetic 
field. Turbulence was shown to be able to prevent 
gravitational collapse (cf Chandrasekhar 1958; Arny 
1971; Bonazzola et al. 1987, 1992) in the 2-D numer- 
ical simulations by Leorat, Passot & Pouquet (1990). 
Vazquez-Semadeni, Passot & Pouquet (1995) mod- 
eled the galactic disc, on the scale of 1 Kpc, as a 
turbulent self-gravitating flow. They simulate a two 
dimensional turbulent flow that is forced by the en- 
ergy released by star formation (expansions of HII re- 
gions), and found that the main mechanism of cloud 
formation is the turbulent ram pressure, rather than 
gravity. They were not able to form self gravitating 
clouds, due to limitations in the thermal modeling, 
and the consequent low density contrast. 

A three dimensional description of a magnetized 
self-gravitating cloud was given by Carlbcrg & Pu- 
dritz (1990), and was used to simulate molecular emis- 
sion spectra by Stenholm & Pudritz (1993). Carlberg 
& Pudritz found that the magnetic field and hydro- 
magnetic waves can support the cloud against grav- 
ity. The clouds contract, because of ambipolar dif- 
fusion, on a time scale of approximately four free- 
fall times. These simulations do not solve the MHD 
equations, but instead make use of a 'sticky particles' 
code. Energy is injected in the form of a spectrum 
of Alfven waves, and the outcome of the computation 
is dependent on the spectral index, that is a free pa- 
rameter. This way of forcing the particles is rather 
unphysical, because an arbitrary spectrum of Alfven 
waves is imposed, instead of being obtained as a re- 
sult of the simulated magneto- hydrodynamics. Pas- 
sot, Vazquez-Semadeni & Pouquet (1995) introduced 
the magnetic field in their previous two dimensional 
model for the galactic disc (Vazquez-Semadeni, Pas- 
sot & Pouquet 1995), and obtained a flow with rough 
cquipartition of kinetic and magnetic energy, proba- 
bly in rough cquipartition also with the mean thermal 
energy. The same simulation, and others with larger 
density contrast and resolution, have been studied by 
Vazquez-Semadeni, Ballesteros-Paredes & Rodriguez 
(1997), who were able to reproduce the observed rela- 
tion between line- width and size (Larson 1981). Gam- 
mie and Ostriker (1996) solved the MHD equations in 
a slab geometry, including self-gravity By forcing the 
flow with a nonlinear spectrum of MHD waves, they 
were able to prevent the gravitational collapse of their 
one dimensional cloud model. 



Apart from the intentional exclusion of gravity, we 
differ from previous studies of turbulent flows in sev- 
eral respects. First of all we have solved the MHD 
equations in three dimensions, while all previous so- 
lutions are in two or one dimensions. In MHD the 
dimensionality of the flow has a fundamental impor- 
tance, because it determines the topological freedom 
of the magnetic field. A price that has to be paid 
for using fully three-dimensional simulations is a sig- 
nificant reduction of the number of points in each 
spatial direction. Thus, while our 128 3 runs corre- 
spond to ~ 2 10 6 degrees of freedom, two dimensional 
simulation such as those of Vazquez-Semadeni et al. 
(1997), that use up to 800 points in each direction, 
still need less than a third of the computing resources 
of our three dimensional simulations. Studies of vary- 
ing dimensionality are thus complementary, and may 
contribute to uncover different aspects of the system 
under study. 

An important property of our models is the high 
Mach numbers of the flows (rms up to ~ 30, maximum 
~ 60), that allows us to investigate the super-sonic 
regime of interstellar medium flows on rather large 
scales, up to about 40 pc. Previous works have nor- 
mally been limited to smaller rms Mach number. For 
example some values of the initial rms Mach number 
in decaying flows are: Mi — 4 in the two dimensional 
experiment by Passot, Pouquet & Woodward (1988); 
Aii = 1 in three dimensional runs by Porter, Pou- 
quet & Woodward 1992; Mi = 10 in one dimensional 
simulations by Gammie & Ostriker (1996); Mi = 5 
in MacLow et al. (1998) in three dimensions. Other 
values of rms Mach number in driven numerical ex- 
periments are: M = 3 in three dimensional simu- 
lations by Lee, Lele & Moin (1991); M = 3 in nu- 
merical experiments by Scalo et al. (1998); M = 5 
in three dimensional simulations by Stone, Ostriker 
& Gammie (1998). The two dimensional simulations 
by Vazquez-Semadeni, Passot & Pouquet (1996) have 
large Mach numbers because of the cooling they in- 
clude. 

While previous models of magnetized clouds fo- 
cused on the role of the magnetic field as opposed 
to gravity, and therefore have assumed a magnetic 
pressure much larger than the gas pressure, the main 
purpose of the present work is to show that MCs are 
well described as flows with lower magnetic pressure 
than previously assumed. 

In this work we report the results of numerical ex- 
periments in two opposite regimes; Ma 3> 1 and 
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M.a ~ 1, where M.a is the initial Alfvenic Mach 
number -the ratio of the rms velocity to the Alfven 
speed. We argue that the observations of magnetic 
field strengths are consistent with a scenario where 
the mean magnetic pressure is dynamically low; i.e., 
where M\ 3> 1. Super-sonic and super-Alfvenic 
magneto-hydrodynamic (MHD) flows develop a very 
intermittent spatial distribution of the magnetic en- 
ergy, such that when the field is detected at a favor- 
able position, its estimated strength is far higher than 
the mean field strength. Therefore, dense cores with 
sub-Alfvenic velocity dispersion can be generated in 
super-Alfvenic flows, in agreement with the observa- 
tions. 

2. The experiments 



The stilidy ot the dynamics ot iVlCs belongs to the 



field of random and super-some MriD flows. The 
Reynolds number and the magnetic Reynolds num- 
ber in MC flows are very large. The random nature 
is therefore a basic feature of the dynamics of these 
flows, and requires an appropriate description. 

A realistic description of the dynamics of molecu- 
lar clouds have to be based on the numerical solution 
of the compressible MHD equations in three dimen- 
sions, in a regime of random and highly supersonic 
motions. We solve the compressible MHD equations 
in such a regime, at a numerical resolution of 128 3 , 
with a code designed for turbulence and MHD tur- 
bulence experiments (Nordlund, Galsgaard & Stein 
1994; Stein, Galsgaard & Nordlund 1994; Galsgaard 
& Nordlund 1996; Nordlund, Stein & Galsgaard 1996; 
Nordlund & Galsgaard 1997), specifically adapted to 
be able to deal with very strong shocks and very large 
density contrasts. 

Although we have already developed a version of 
the code with the inclusion of the gravitational force, 
all the experiments were run without gravity, for the 
following reasons: 

• We are interested in studying if gas dynamics 
and magneto-hydrodynamics alone can shape 
distributions that are similar to those observed, 
even without invoking gravity. 

• The observed motions have velocities compara- 
ble with the virial velocity, or larger, on a range 
of scales, and the clouds are not free-falling. 

• If the results of our experiments are discussed 
only up to a time shorter than or comparable 



to the dynamical (or free-fall) time, all our con- 
clusions remain basically unchanged. This time 
is about a few million years on a scale of 10 pc, 
and clouds are not supposed to live much longer 
than that, before star formation takes place and 
becomes energetically important. 

We note that the fact that clouds are observed to 
have velocities comparable to virial velocities, and 
are observed to not be free- falling, are not accept- 
able justifications for neglecting self-gravity; a com- 
plete model of interstellar clouds should instead ex- 
plain why they are not free-falling even in the pres- 
ence of gravity. However, by concentrating first, in 
the present paper, on a situation where gravity is ne- 
glected, we provide a point of reference for subsequent 
studies, where we do include self-gravity (Padoan 



et al., 1999b|) 



2.1. The equations 

We solve the compressible MHD equations: 
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plus numerical diffusion terms, and with periodic 
boundary conditions, v is the velocity, B the mag- 
netic field, f an external force (= in these particu- 
lar experiments), and P — pT is the pressure at T « 
const. 

Conditions in the cold molecular clouds that we 
are modeling are such that an isothermal approxima- 
tion is adequate; the radiative heat exchange is so 
efficient that the temperature remains low in most 
places. Even if the temperature momentarily in- 
creases in shocks, the subsequent cooling is rapid, and 
the result is shock structures that are qualitatively 
and quantitatively similar to isothermal shocks. 
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We have thus used isothermal conditions in most 
of our runs and have verified that this is appropri- 
ate, by rerunning segments of some experiments us- 
ing the full energy equation. No significant change 
of the statistics was found and, since using the full 
energy equation increases the cost of the experiments 
considerably (the strong cooling required to maintain 
a low temperature forces a much smaller time step), 
we performed most of the experiments at constant 
temperature. 

The absence of an explicit resistivity r\ in the in- 
duction equation corresponds to an assumption of flux 
freezing on well resolved scales. The code uses shock 
and current sheet capturing techniques to ensure that 
magnetic and viscous dissipation at the smallest re- 
solved scales provide the necessary dissipation paths 
for magnetic and kinetic energy. As shown by Gals- 
gaard & Nordlund (1996, 1997), dissipation of mag- 
netic energy in highly turbulent, compressible MHD 
plasmas occurs at a rate that is independent of the 
details of the small scale dissipation. In ordinary hy- 
drodynamic turbulence the corresponding property is 
one of the cornerstones of Kolmogorov (1941) scaling. 

The works by Galsgaard & Nordlund (1996, 1997) 
refer to compressible flows, but since the ratio of mean 
kinetic energy to mean magnetic energy is small, mo- 
tions arc in practice roughly incompressible, at least 
in the directions perpendicular to the local magnetic 
field. In the case referred to in the present paper 
as Model A, the mean kinetic energy is instead large 
compared to the mean magnetic energy, and the sit- 
uation might be rather different, and more similar to 
Burger's turbulence than both classical Kolmogorov 
turbulence, and the situation studied by Galsgaard 
and Nordlund. However, Burger's turbulence is also 
characterized by having a dissipation that is, in the 
limit, independent of the value of the viscosity (cf. 
Saffman 1968, p. 485, sec. 6). 

We have not included ambipolar diffusion in any of 
the driven experiments. The value of the rms Mach 
number in the driven experiments is comparable to 
the value estimated in molecular clouds on a scale 
range from 3 pc to 50 pc. On such scales, ambipolar 
diffusion occurs in a time significantly longer than 
the dynamical time, as recently shown by Myers & 
Khcrsonsky (1995) and is expected to be of secondary 
importance (see also Elmegreen & Fiebig 1993). 

On sufficiently small scales and high densities, am- 
bipolar diffusion can have significant effects (Zweibel 
& Brandenburg 1994, 1997; Tagger et all995). Since 



the decay runs have lower rms Mach number (smaller 
scale), and since ambipolar diffusion is considered 
an important dissipation mechanism in magnetized 
clouds (McKee et al. 1993), the ambipolar drift has 
been included in some of the decaying runs. 

2.2. The code 

The code solves the compressible MHD equations 
on a 3D staggered mesh, with volume centered mass 
density and thermal energy, face centered velocity and 
magnetic field components, and edge centered electric 
currents and electric fields (Nordlund, Stein & Gals- 
gaard 1996). 

The original code works with "per-unit-volumc" 
variables; mass density, momenta, and thermal energy 
per unit volume. In the super-sonic regime relevant 
in the present application, we found it advantageous 
to rewrite the code in terms of "per-unit-mass" vari- 
ables; lnp, u, and E = With these variables, the 
time evolution of all variables is governed by equa- 
tions of the type 



D f df 

— = — + v • Vf = 

Dt at J 



(6) 



i.e., equations that specify the time rate of change fol- 
lowing the motion. These are better conditioned than 
the divergence type equations that result from us- 
ing per-unit-volume variables (the large — order M 2 — 
density variation in isothermal shocks cause the per- 
unit-volume fluxes to vary over several orders of mag- 
nitude). 

We use spatial derivatives accurate to 6th order, 
interpolation accurate to 5th order, and Hyman's 3rd 
order time stepping method (Hyman 1979). 

In order to minimize the viscous and resistive in- 
fluence on well resolved scales, we use monotonic 3rd 
order hyper-diffusive fluxes instead of normal diffu- 
sive fluxes, and in order to capture hydrodynamic 
and magneto-hydrodynamic shocks we add diffusiv- 
ities proportional to the negative part of the velocity 
divergence, and resistivity proportional to the nega- 
tive part of the cross-field (two-dimensional) velocity 
divergence. Further details of the numerical methods 
are given by Nordlund, Stein & Galsgaard (1996) and 
Nordlund & Galsgaard (1997). 

2.3. Weak and strong magnetic field 

For the purpose of this work we have run several ex- 
periments, that can be divided into two main groups: 
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i) super-Alfvenic runs (models A), with Mas ~ 10 
(in one case Mas ~ 30), and 0i > 1; ii) equipar- 
tition (of kinetic and magnetic energies) runs (mod- 
els B), with Mas ~ 1 and <C 1; where Mas is 
the initial value of the Alfvenic Mach number (the 
rms velocity of the flow divided by the Alfven ve- 
locity) and Pi = (P g /P m )i is the initial ratio of gas 
pressure to magnetic pressure. Notice that the mag- 
netic energy grows considerably in the beginning of 
the super-Alfvenic runs, and therefore the value of 
Ma can quickly decrease by a factor of two and the 
value of j3 by a factor of four. 

Both models A and B are divided in two sub- 
groups: decaying runs (simply A and B), and ran- 
domly driven runs (Ad and Bd) . Models of the same 
type can also differ by the initial rms Mach number of 
the flow (the rms velocity of the flow divided by the 
speed of sound), or by the inclusion or not of ambipo- 
lar diffusion. A list of the models is given in Table 1, 
together with the initial values of their parameters. 

In all experiments the initial density is uniform, 
and the initial velocity is random. We generate the 
velocity field in Fourier space, and we give power, 
with a normal distribution, only to the Fourier com- 
ponents in the shell of wave- numbers 1 < kL/2w < 2. 
We perform a Helmholtz decomposition, and use only 
the solenoidal component of the initial velocity. The 
initial magnetic field is uniform, and is oriented par- 
allel to the z axis: B = Bqz. 

3. Results 

In this section some properties of the models that 
may be related to observational quantities are dis- 
cussed. The initial values of some physical parame- 
ters of these models are listed in Table 2. There are 
only two numerical parameters in the models: M and 
Ma- In order to rescale the models to physical units, 
we use the following empirical Larson's relations (Lar- 
son 1981): 
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where a temperature T ~ 10 K is assumed, and 



L 



(n) = 2.0 x 10-M 

lpc 



(7) 



(8) 



the value of L, the linear size of the model cloud, 
and of (n) , the average number density. The physical 
unit of velocity in the code is the isothermal speed of 
sound, C s , and the physical unit of the magnetic field 
is C s (47r(p})5 ( C gs). 

3.1. Dissipation of supersonic motions 

Molecular cloud lifetimes are estimated to be of 
the order of a few dynamical times, while super- 
sonic motions should dissipate their kinetic energy 
in shocks in about one dynamical time. One of the 
original motivations for models of magnetized clouds 
with equipartition of kinetic, magnetic and gravita- 
tional energies was the belief that MHD waves decay 
on a longer time-scale than random super-sonic and 
super-Alfvenic flows, and therefore could offer an ex- 
planation for the long lifetime of molecular clouds. 
Zweibel & Josafatsson (1983) and Elmegreen (1985), 
noticed that decay of MHD waves can be rather short, 
but their studies are based on the ambipolar diffusion 
dissipation. Ambipolar diffusion is in fact believed 
to be the most important dissipation mechanism in 



that is equivalent to a constant mean surface density. 
With these two relations, the value of M determines 



magnetized clouds (McKce et al., 1993). 

We have run decaying experiments both with and 
without ambipolar diffusion, from which we can see 
that the inclusion of ambipolar diffusion does not 
change very much the decay time-scale. In our de- 
caying equipartition models, the most important en- 
ergy dissipation occurs because motions along the 
magnetic field lines produce strong shocks, just as 
in super-Alfvenic models. This was found also by 
Gammie and Ostriker (1996). The ratio of the decay 
time to the instantaneous dynamical time is plotted 
in Fig. 1 for models without ambipolar diffusion, and 
in Fig. 2 for similar models, but with the inclusion of 
ambipolar diffusion. Decay times of kinetic, magnetic 
and total energies are plotted versus the time in units 
of the initial dynamical time, Td yn = L/a v . After a 
short initial transient (of the order of 0.5Td yn because 
the initial velocity field has power up to kL/2-K = 2), 
the decay time is in all cases comparable to the instan- 
taneous dynamical time. Thus, the strength of the 
magnetic field, and also the importance of ambipolar 
diffusion, do not significantly affect the time-scale of 
decay of the flow. 

A simple interpretation of the insensitivity of the 
level of dissipation to details of the dissipation mech- 
anism is that it is related to the origin of the scal- 
ing behavior of turbulent flows. We observe a similar 
behavior in our super-sonic experiments as in experi- 
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ments with sub-sonic turbulence; a power law energy 
power spectrum is established already after a fraction 
of a dynamical time. The scaling is such that the 
decay of energy does not depend on the particular 
dissipation process that occurs on the small viscous 
and resistive scales, but is rather determined by the 
large scale kinematics, characterized by the dynami- 
cal time-scale. 

It may be concluded that equipartition models do 
not offer any advantage as far as the energy decay 
is concerned, compared with super-Alfvenic models. 
Mac Low (1998)| and |Stone et al. (1998)| have recently 



shown that the decay time in equipartition models can 
be even shorter than in super-Alfvenic models. The 
problem of the lifetime of molecular clouds is not to be 
solved by the intervention of a strong magnetic field, 
but more likely by continuous driving of the turbu- 
lence. Cascading energy from distant supernova ac- 



ti vity is a likely source of galactic turbulence ( Korpi 



al, 1998b| ; [Korpi ct al., 1998a| ; |Korpi et al., 1999| ) 



et 

Also, bipolar outflows from stars formed inside molec- 
ular clouds can provide significant amounts of energy 
(Reipurth, Bally & Devine, 1997). Moreover, cloud 
life times are not very long (Blitz & Shu 1980), and 
could perhaps even be comparable with their dynam- 
ical time (Elmegreen & Efremov, 1998). 

In this connection it should be pointed out that, 
even though energy sources on the average replenish 
turbulence in clouds, the driving may be intermittent 
in time as well as in space, and turbulence in individ- 
ual clouds may thus, at times, be decaying. 

3.2. Probability Distribution of Magnetic Field 
Strength 

The main difference between the super-Alfvenic 
model and the equipartition model is the fact that in 
the latter the kinetic energy is not sufficiently large 
to compress the magnetic field significantly, while in 
the former very strong density enhancements pro- 
duce large local enhancement of the magnetic field 
strength. Moreover, expansions generate large voids 
with magnetic field strength much smaller than the 
mean field in the super-Alfvenic model, but not in 
the equipartition model. 

As a result, runs of type A produce very inter- 
mittent probability distributions of the field strength, 
with a roughly exponential tail. Models of type B, 
instead, are characterized by approximately Gaus- 
sian probability distributions of the field strength. In 



Fig. ||, the field strength distributions from the last 
snapshots of model Ad2 (thick line) and model Bdl 
(thin line) are plotted. In model Ad2, (B) = 4.5 ^G, 
but values of field strength up to 100 fiG are reached 
in high density regions. The values of B span almost 
four orders of magnitude, 0.03 ^G< B <100 ^G, in 
model Ad2, while they span only one order of magni- 
tude in model Bdl, 4.5 fiG< B <46 ^G. 

In the super-Alfvenic model most of the mass con- 
centrates in regions of large gas density (see below), 
where the field strength is also large. Because of this, 
approximately 0.5% of the total mass of the system 
contains a field 10 times stronger than the mean value. 
For example, the last snapshot of model Ad2 corre- 
sponds to molecular clouds on the scale of 15 pc, with 
(n) = 133 H2Cm~ 3 . It may be taken, therefore, as a 
model for a molecular cloud of about 2 x 10 4 M . 
Since (B) = 4.5 /jG, one clump of about 100 M can 
be found in the cloud, with B ~ 45 /jG, according to 
model Ad2. 

The examples show that the formation of cores 
with rather large field strength is predicted in the 
super-Alfvenic model, even if the mean magnetic field 
strength in the cloud is quite small, and comparable 
with the mean galactic field strength. As a result, 
Zeeman splitting observations should not easily detect 
the magnetic field in clouds, apart from inside par- 
ticularly favorable high density clumps. Conversely, 
if the magnetic field is detected in a cloud core, it is 
likely that the mean field in the cloud is much smaller 
than the estimated field strength in that core, if the 
cloud dynamics is described by a random super-sonic 
and super-Alfvenic flow, as in model Ad2. 

3.3. Probability Distribution of Gas Density: 
Stellar Extinction 

In a highly radiative gas, the most important ef- 
fect of supersonic random flows is the fragmentation 
of the medium into filaments and clumps with very 
large density contrast. Supersonic turbulence can in- 
deed be the reason for the observed fragmentation of 
molecular clouds (eg Padoan & Nordlund 1998). It is 
interesting therefore to quantify the statistics of the 
density distribution that arise from supersonic tur- 
bulence. Vazquez-Semadeni (1994) found that the 
density distribution in his two dimensional numeri- 
cal simulations was consistent with a Log-Normal. 
Padoan, Nordlund & Jones (1997) also found density 
distributions consistent with a Log-Normal in their 
three dimensional numerical simulations, while Scalo 
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et al. (1998) studied the effect of a polytropic equa- 
tion of state on the density distribution. Nordlund 
& Padoan (1998) used a new set of numerical exper- 
iments to show that even polytropic flows produce a 
density distribution similar to the Log-Normal, but 
slightly skewed. They also argue that such a distri- 
bution is due to the fact that the dynamics does not 
depend on the mean density. 

Observationally it is difficult to measure the gas 
volume density in molecular clouds, but stellar ex- 
tinction measurements can provide information about 
the projected gas density in clouds. Lada et al. 
(1994) suggested the use of stellar extinction measure- 
ments as a test for models of the structure of molec- 
ular clouds. Padoan, Jones, & Nordlund (1997) have 
shown that near-infrared stellar extinction measure- 
ments can be used to infer the three dimensional prob- 
ability distribution of the gas density in dark clouds. 
They have shown that there is qualitative and quan- 
titative agreement between the inferred properties of 
the three dimensional density distribution in the dark 
cloud IC5146 (Lada et al. 1994), and the properties 
of the three dimensional density distribution in their 
experiments of random supersonic flows. 

The comparison between numerical models and ob- 
servational data is based on the plot of the dispersion 
of the extinction measurements in cells, versus the 
mean extinction in the same cells (Lada et al. 1994). 
In Padoan, Jones & Nordlund (1997), the theoreti- 
cal plots are calculated starting from random density 
fields of given statistics and power spectra. Here we 
generate the same plot, but calculated directly from 
the density field of the experiments. A random dis- 
tribution of stars is generated and it is assumed that 
most stars are in the background (that is the model 
cloud is between the stars and the observer) . The vi- 
sual extinction is proportional to the column density 
of the cloud sampled by the line-of-sight to any given 
star: 



Ay I mag 



N(H + H 2 



(9) 



2 x 10 21 cm- 2 
(Bohlin et al. 1978) . A uniform grid is used to sub- 
divide the whole area into cells, each containing an 
average number of five stars, as in Lada et al. (1994). 
The mean extinction Ay and its dispersion a(Ay) 
are measured inside each cell, and cr(Ay) is plotted 
against Ay. 

We have performed this test with all our models. 
In Fig. H we show the results for the models Ad2 



and Bdl. These two models both have very large 
rms Mach number and differ from each other because 
model Ad2 is super-Alfvenic, while model Bdl has 
approximate equipartition of magnetic and kinetic en- 
ergy. We use the last snapshots from the two runs, 
where the values of Mach number and Alfvenic Mach 
number are M ss 15 and Ma ~ 5 respectively, for 
model Ad2, and M « 12 and Ma ~ 1 for model 
Bdl. According to the Larson's relations we use in 
the present work to rescale our numerical models to 
physical units (Eqs. (7) and (8)), these snapshots 
correspond to molecular clouds on the scale of 14 pc 
(Ad2) and 9 pc (Bdl). Using the data by Lada et al. 
(1994) and by Dobashi et al. (1992), one may infer an 
rms Mach number M 10 for the dark cloud IC5146, 
comparable to the Mach number values in the simula- 
tions. However, the particular scale is not important 
as for the mean or total column density of model or 
real clouds. This is because the Larson's relations are 
such that the column density of clouds is independent 
of size. The scale is instead important for determin- 
ing the slope of the upper envelope of the Ay-a(Ay) 
plot. As shown in Padoan, Jones & Nordlund (1997), 
that slope is roughly proportional to the standard de- 
viation of the three dimensional density distribution, 
which is also proportional to the Mach number of the 
flow. So the slope is expected to be larger on larger 
scales (higher rms Mach number). 

It is clear for Fig. |^ that the super-Alfvenic model 
reproduces very well the observed plot, while the 
equipartition model (Bd4) does not. Although model 
Bd4 has an rms Mach number comparable to the one 
in the cloud IC5146, the values of Ay and a(Ay) in 
the model span a much smaller range than in the ob- 
servations. This is also illustrated by the histogram of 
extinction (column density), shown on the left panels 
of Fig. H The equipartition model behaves like an 
elastic medium in the directions perpendicular to the 
magnetic field lines, and this reduces very much its 
ability to generate a sufficiently complex density field. 
In particular, while compressions on all three dimen- 
sions can generate a topology with large holes and 
adjacent lines-of-sights with very large fluctuations 
in extinction (model Ad2), compressions in only one 
dimension (model Bdl) can do that only to a smaller 
extent. 
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3.4. Synthetic spectra: line width versus in- 
tegrated temperature 

Most of the data from dark clouds consist of molec- 
ular emission line observations. The intensity of a 
particular emission line depends on many factors such 
as gas density, kinetic temperature, radiation field, 
and chemical abundances. Moreover, the line pro- 
files are strongly dependent on the gas kinematics, 
but in a rather complex way, because only radial ve- 
locities are available, and both line widths and in- 
tensities are quantities projected along the line-of- 
sight. It is very difficult to extract reliable informa- 
tion about the three dimensional structure and kine- 
matics of molecular clouds directly from observational 
data, and making use of simple models assuming lo- 
cal thermodynamic equilibrium (LTE), smooth ve- 
locity fields, and oversimplified density distribution. 
Part of the problem is due to the complexity of the 
cloud structure and of the internal random motions 
in clouds. 



Padoan et al. (1998b) have recently calculated 
maps of synthetic molecular spectra that may be used 
to infer intrinsic properties of molecular clouds by 
comparison with observed spectra. Details of the 
method can be found in Padoan ct al. (1998b)| . Here 
we simply remark that the synthetic spectra are ob- 
tained with a non-LTE Monte Carlo radiative transfer 
code (Juvela 1997), and the radiative transfer calcu- 
lations are performed on a cloud model that is sim- 
ilar to the ones used in present work. Images of 
synthetic molecular map s of different transitions and 
molecules are shown in Padoan et al. (1998b) , to- 



gether with statistics of the spectral line profiles. In 
that work, a super-Alfvenic flow was used, and was 
shown to reproduce observed properties of clouds (see 
also Padoan et al., 1998a ) . In the present work we are 
interested also in the comparison between the super- 
Alfvenic and equipartition models, and we have there- 
fore solved the radiative transfer problem using cloud 
models Ad2 and Bd4 (the radiative transfer calcula- 
tions have been kindly provided by Mika Juvela). No- 
tice that although the initial Mach number of model 
Ad2 is almost 30, it later decreases (due to its initial 
conditions and not strong enough random driving) , to 
a value comparable to the one in model Bd4. We use 
the last snapshots from each model, which correspond 
to the physical regime in typical molecular clouds on 
the scale of about 10-15 pc. 

Comparisons of molecular spectra between the two 



models have shown that while the two cases are 
rather similar in many of their observable proper- 
ties, they are very different in the relation between 
J=1^0 13 CO line width and J=1^0 13 CO integrated 
antenna temperature. The integrated antenna tem- 
perature is roughly proportional to the gas column 
density (cf. Dickman 1978, and Padoan et al., 1999a, 
for a critical discussion based on non-LTE radiative 
transfer), and its histogram is shown in Fig. |^. The 
equipartition model has a narrower distribution than 
the super-Alfvenic model, consistent with the dis- 
cussion of dust extinction in the previous subsection. 
The histogram of line width is also narrower for model 
Bd4 that for model Ad2 (not shown), although the 
spectra averaged over all synthetic maps are almost 
identical (same line width) . The relation between the 
line width and the integrated antenna temperature is 
also quite different in the two cases. The equiparti- 
tion model shows almost no growth of line width with 
increasing J=1^0 13 CO antenna temperature, while 
in model Ad2 the growth is very significant, as shown 
in Fig.| 

Heyer, Carpenter & Ladd (1996) find that in the 
star forming giant molecular clouds they observed the 
J=1^0 12 CO line width grows with J=l->0 13 CO in- 
tegrated antenna temperature (roughly proportional 
to column density), and discuss the possibility of us- 
ing this property to infer the dynamic importance of 
the magnetic field in the cloud motions. They ar- 
gue that, if the motions observed in molecular clouds 
were Alfven waves, like in theoretical models of dark 
clouds (for example in Carlberg & Pudritz 1990), then 
the line width should decrease with column density. 
Based on this, they conclude that the motions in the 
clouds they observed are not consistent with Alfven 
waves. We find that the J=1^0 13 CO line width 
in the equipartition model is in general almost inde- 
pendent of column density (J=l-^0 13 CO integrated 
antenna temperature), and in some cases even de- 
creases with column density, as predicted by Heyer, 
Carpenter and Ladd (1996) for the J=1^0 12 CO line 
width. The equipartition model is therefore in conflict 
with the observations, while the super-Alfvenic model 
shows a relation between line width and column den- 
sity that is similar to the observational result. 

In order to better constrain the importance of the 
magnetic field in the dynamics of molecular cloud in- 
ternal motions, it would be useful to measure line 
width versus integrated temperature in many differ- 
ent clouds. We expect that the magnetic field could 
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play an important role in some cores with large field 
strength, and this could appear in the relation be- 
tween J=l— >0 13 CO line width and column density. 
Cores with strong magnetic field (predicted in the 
present super- Alfvenic scenario for molecular cloud 
dynamics) have already been identified by OH Zee- 
man splitting measurements, and should be the target 
of such a study. 

3.5. The B-n relation 

As discussed above, the equipartition model gen- 
erates a smaller density contrast than the super- 
Alfvcnic model, and also generates much smaller fluc- 
tuations in the magnetic field strength. The relation 
between field strength and gas density (B-n relation) 
is therefore expected to be rather different in the two 
models. The B-n relation in the models may be com- 
pared with observational results. 

In regions of maser emission, at densities of about 
n = 10 7 cto~ 3 , a field strength of the order of B = 
10 3 — 10 4 ^G is observed, while in regions of molecu- 
lar emission, with approximately n — 10 2 — 10 3 cm -3 , 
the field is found to be of the order of B = 10/iG 
(eg Myers & Goodman 1988). A relation of the type 
B oc n°- 3 ~ - 6 may be deduced from the observations 
(eg Troland, Crutcher & Kazes 1986; Heiles 1987; Du- 
dorov 1991), but it is quite uncertain, especially in the 
light of the above discussion about the intermittency 
of the distribution of the magnetic energy. 

We have updated the B-n relation by Troland & 
Heiles (1986) with more recent measurements: the HI 
Zeeman observations by Verschuur (1995); the OH 
and CN Zeeman splitting measurements by Crutcher 
et al. (1993, 1996, 1998); the OH maser measure- 
ments by Johnston et al. (1989); the H2O maser 
measurements by Fiebig & Gusten (1989). 

In Fig. 0, observational detections of the field- 
strength and upper limits are plotted versus the es- 
timated value of the gas density. Although most of 
the attempts to detect the Zeeman effect result in 
non-detections, and therefore in upper limits to the 
field-strength, there are some detections of B also in 
the range of density found in molecular clouds. B is 
of course detected preferentially in the regions where 
it is particularly strong, which are hardly representa- 
tive of the typical physical conditions inside molecular 
clouds. Fig. ^ shows a clear power-law relation be- 
tween B and n, defined by the regions with strongest 
magnetic field. At the same time, there are many es- 



timated upper limits (over 100 listed in publications 
and only summarized in Fig. ^) that very likely probe 
many regions with values of B weaker than the ones 
detected in regions of similar densities. The observa- 
tions therefore clearly indicate that the B — n relation 
has a large scatter, although a power-law upper en- 
velope is well defined over approximately 10 orders of 
magnitude in gas density and 5 orders of magnitude 
in field-strength. The upper envelope is roughly 

^^(h^) ■ < 10 > 

and the scatter at n = 10 3 cm~ 3 is at least two orders 
of magnitude, in the approximate range 1 — 100 /iG. 

A super- Alfvenic flow (eg model Ad2) naturally de- 
velops a B — n relation with a power law upper enve- 
lope and a large scatter, and also provides a large den- 
sity contrast at the same time. The B-n relation is 
initially almost linear (see below for an explanation), 
and flattens with time. The upper envelope of the 
relation also flattens with time, and after about one 
dynamical time it stabilizes to B oc n 0A , as in the ob- 
servations. In Fig. ^ contour lines obtained from the 
B-n scatter plot in model Ad2 (thick lines) are shown, 
together with the observational data. Model Ad2 cor- 
responds to typical molecular clouds on the scale of 
20 pc (rms Mach number « 20). The particular snap- 
shot used in Fig. [?] has M = 15 and (B) = 4.5 uG. 
Although the mean field strength is so low (of the 
order of the mean galactic field), the fluctuations of 
field strength are every large (0.03 < B < 100 /LtG), 
and the largest values of the magnetic field strength 
(upper envelope) at each density are very close to the 
observational values. In model Ad2 there is a clear 
power-law upper envelope over about 4 orders of mag- 
nitude in density and almost 2 in B, which matches 
the observations nicely. The scatter in the B-n rela- 
tion of model Ad2 is also consistent with the observa- 
tions. In order to achieve the densest values of field 
strength estimated in molecular clouds using CN Zee- 
man splitting and masers, higher numerical resolution 
would be necessary, or perhaps self-gravity should be 
included in the calculations. 

The result for an equipartition model (Bdl) is also 
overplotted in Fig. ^ (thin lines). The model is in- 
consistent with the observations in two ways: i) the 
ranges of density and field strength values are too 
small (see the two subsections above) and ii) neither 
the upper envelope nor the mean field strength at each 
density define a B-n relation similar to the observa- 
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tional one. In particular, most observational upper 
limits on the field strength, and a few detections, are 
completely inconsistent with the equipartition model. 

While super sonic and super- Alfvenic turbulence 
naturally explains the results of Zeeman splitting 
measurements, including the B-n relation and its 
scatter, such results are not easily explained by the 
alternative equipartition model. 

The lack of correlation between B and n in model 
Bdl is due to the fact that most of the density en- 
hancement occurs by convergent flows along the field 
lines, and these motions do not affect the magnetic 
field. 

One could possibly argue that the lack of corre- 
lation of B and n in model Bdl is consistent with 
the majority of observations, because often no field 
is detected, and hence nothing can be said about the 
correlation. However, as further discussed in the next 
Section, the upper limits of the non-detections speak 
against this interpretation. Also, the cases where a 
field is detected would then have to be explained with 
ad hoc arguments, rather than as a natural part of a 
statistical distribution. 

3.6. Cloud and flow topology 

An understanding of the spatial structure of the 
density field in dark clouds is very important for a 
correct interpretation of observational data, and for 
the formulation of a number of physical models. In 
this subsection we provide our interpretation of the 
dynamical phenomena occurring in the simulations, 
based on extensive three dimensional browsing of the 
data, and illustrated here by images of the density 
field and plots of statistical correlations of the velocity 
and magnetic fields. 

The snapshots in Fig. || and Fig. ^| give an idea 
of the dimensionality of the structures in the den- 
sity field. The topology of the density field in ex- 
periment Ad2 (Fig. |^) has a clear evolution in time. 
In the very beginning, until t m 0.6td yn ,o, the den- 
sity grows predominantly in sheets. These are the 
fronts of blobs of coherent motion, advancing at su- 
personic velocity. Later, these fronts start to intersect 
each other, and the density increases especially in fil- 
aments (at the intersections of fronts). The evolu- 
tion continues with the intersection of filaments into 
knots, at t ~ 1.5i<2yn,o- The fully developed topol- 
ogy, at t = 2.0td y m is characterized by both filaments 
and knots (cores). Also in experiment Bdl (Fig. H) 



the density grows initially in sheets, but the transi- 
tion from sheets to filaments is mainly due to motions 
along field lines, rather than to intersection of sheets. 
In the equipartition model Bdl, in fact, density en- 
hancements cannot move freely in the three space di- 
rections, but only along magnetic field lines. It is 
therefore typical that density enhancements such as 
sheets are torn apart, by motions along field lines, 
into filaments approximately aligned with the mag- 
netic field. The approximate alignment of the fil- 
aments is visible in the left panel of Fig. |[ Later 
in time, the density structure becomes more complex 
(right panel of Fig. ||), but still shows some evidence 
of the alignment of the density filaments. 

The evolution of the magnitudes of the mass den- 
sity and the magnetic field may be discussed with ref- 
erence to Lagrangian version of the continuity equa- 
tion, 

^-V-t, (U) 
and the scalar induction equation 



Dln|S| 
Dt 



-V- 



B ■ (B ■ V)vj 



(12) 



where v_l- is the velocity perpendicular to the mag- 
netic field, and —V • vi usually is the dominant term 
on the RHS. 

Although — V • v vanishes for the solenoidal initial 
condition, the supersonic motions rapidly lead to the 
formation of shock fronts, where the local value of 
— V • v is large and positive because of the disconti- 
nuity in the velocity perpendicular to the shock. 

The initially homogeneous magnetic field is carried 
along by the perpendicular components of the veloc- 
ity field, and is hence also collected into sheets, except 
at those rare locations where the initial field happens 
to be strictly parallel to the velocity field. This ex- 
plains why sheets initially form in both mass density 
and magnetic flux density, and why the B-n relation 
initially has an exponent close to unity. 

Note that the usual argument for a slope of 2/3, 
that applies to isotropic and non-shocking compress- 
ible motion does not apply here, because of the devel- 
opment of discontinuities. In term of Eqs. |ll| and [l^, 
the 2/3 follows if — V • vj_ typically picks up two of 
three statistically equivalent contributions to —V • v. 
However, at a shock, the divergence is dominated by 
the derivative in one particular direction; the one per- 
pendicular to the shock front. The magnetic field that 
is swept into the discontinuity quickly becomes almost 
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parallel to the shock front, because the component in 
the plane of the shock grows exponentially with time. 
Thus, as long as the topology is dominated by sheets, 
the mass density and the magnetic flux grow more 
or less in unison in the sheets, corresponding to an 
exponent in the B — n relation close to unity. 

In the subsequent evolution, there are effects that 
tend to reduce the exponent in the B—n relation. 
First, the non- linear evolution of the initially solenoidal 
velocity field also leads to the development of regions 
of space with a positive divergence, in which both the 
mass density and the magnetic flux density decline. In 
these regions, there is no particular dominance of the 
cross-field divergence, and thus the three-dimensional 
divergence picks up an additional contribution rel- 
ative to the two-dimensional divergence, consistent 
with the classical 2/3 argument outlined above. 

In experiment Bl the motions are mainly Alfven 
waves, and therefore the velocity is predominantly 
perpendicular to the direction of the magnetic field. 
This is illustrated in Fig. [h], where we have plotted 
the histogram of cos(a), where a is the angle between 
v and B. Note that, even though motions across the 
field lines are the most common ones, it is the less fre- 
quent motions along the field lines that dominate the 
dissipation. The motions across field lines are sub- 
ject to magnetic restoring forces, and do not lead to 
substantial density enhancements. 

In experiment Al, the magnetic field is advected 
by the flow, and the stretching of field lines instead 
produces some alignment between v and B, already 
before one dynamical time has passed, as illustrated 
in Fig. 0. 

Alignment between B and v may be caused by two, 
complementary effects: 1) Dynamical alignment is ex- 
pected when the magnetic energy approaches and ex- 
ceeds the kinetic energy; the Lorentz force then forces 
the flow to be predominantly along the magnetic field 
lines. 2) Kinematic alignment occurs when a spatially 
non-uniform velocity field causes stretching of mag- 
netic field lines, and hence a correlation of B and v. 
Pure shear, for example, tends to align an embedded 
magnetic field with the direction of the flow. 

Motions that are aligned with the magnetic field af- 
fect the mass density without affecting the magnetic 
flux density. In particular, the non-linear concentra- 
tion into first sheets and then filaments due to the 
interaction of shock fronts continues into the forma- 
tion of knots in the density field, by the convergence 



of matter flowing along filaments. There is no corre- 
sponding process available to a divergence free vec- 
tor field such as the magnetic field; once the field has 
concentrated into filaments, it cannot concentrate fur- 
ther; the magnetic field in a filament is insensitive to 
flow along the filament. 

In the same way that converging flows along the 
magnetic field may lead to extreme concentrations of 
mass, those regions where the flow is diverging along 
magnetic field lines may lead to extreme rarefactions 
of mass, without affecting the magnetic flux density. 
In B-n scatter plots, this corresponds to the develop- 
ment of more extreme excursions of the mass, relative 
to those of the magnetic flux density and hence a flat- 
tening of the B-n relation with time. 

In model Al, dynamical alignment is at most sig- 
nificant in the few cores that develop a strong (sub- 
Alfvenic) magnetic field in the early evolution of the 
experiment. In scatter plots of B against n most con- 
tributions come from regions where dynamical align- 
ment is unimportant. We thus conclude that the evo- 
lution of the B — n relation in model A, towards a 
smaller exponent with time, is caused by the kine- 
matic alignment of B and v. 

4. Discussion 

It is difficult to get an objective view on the mag- 
netic field strength in dark clouds from the literature. 
The reasons are the following: 

• Negative results from observational programs in 
which detections have been reported are likely 
to remain unpublished. 

• The positions searched for Zeeman splitting 
never represent a statistically meaningful sam- 
ple. Favorable regions are always selected, be- 
cause the observations are very time consuming. 

• The total number of regions in dark clouds, for 
which OH Zeeman observations are published, 
is still small. 

Although all Zeeman splitting measurements are 
biased towards regions where a strong field is expected 
(high density cores) , in many cases the field is not de- 
tected, and upper limits are quite stringent. As an 
example, Crutcher et al. (1993) selected four cores in 
the Taurus dark cloud complex, two in the Libra com- 
plex, two in p Oph, one in the Orion molecular ridge, 
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one position in L889, and the core of Bl (Barnard 1), 
in the Perseus region. The only certain detection is 
in the cloud Bl. For the other regions, the weighted 
average value of the field is +2.7 ± 1.5/iG in Tau- 
rus, -2.1 ± 2.8nG in Libra, +6.8 ± 2.5//G in pOph, 
-0.6 + 2. lfiG in L889, and -4.7+3. 5fiG in L1647. CN 
Zeeman splitting measurements have also provided 
rather low upper limits to the field strength, com- 
pared with estimated equipartition values (Crutcher 
et al. 1996), and also magnetic field detections well 
below equipartition (Crutcher et al. 1998). A rather 
low average field of +9/iG was also found in Cas A by 
Heiles and Stevens (1986). 

Orion A and Orion B are instead two examples 
of detections of intense magnetic field. OH Zeeman 
splitting in Orion A revealed a magnetic field strength 
of B = -125^G (Troland, Crutcher, & Kazes 1986), 
and B = +38//G toward Orion B (Crutcher & Kazes 
1983). 

The OH Zeeman splitting should probe regions 
with n = 10 3 — 10 4 cm" 3 , while the CN Zeeman split- 
ting, should probe regions of n w 10 6 cm" 3 . 

In our model Ad2, the mean field strength is B = 
4.5// G, and the mean density (n) = 1.3 x 10 2 cm -3 ; 
20% of the total mass is in dense clumps, with density 
ten times larger than the mean, n k 3 x 10 3 cm" 3 , 
and field strength five times, or more, larger than the 
mean, B > 25/j, G. Therefore, even if field strengths 
of about 40/iG are detected sometimes in dense cores, 
the mean Alfven velocity in the molecular cloud may 
be just comparable to the sound speed, {v\y/ 2 ~ Cs- 

It may be concluded that model Ad2 is consis- 
tent with the observational estimates of magnetic field 
strength in dark clouds. The particular values of the 
magnetic field used in the comparison here should not 
be taken too literally; the small scale field strengths 
could be larger than estimated by the Zeeman effect if 
the magnetic field is tangled (but notice that in Fig. [?] 
the observational measurements have been multiplied 
by a factor of two, to account for the random orien- 
tation of the magnetic field lines relative to the line- 
of-sight). 

It is difficult to envisage how the measurements 
could be consistent with equipartition models such as 
Bdl or Bd4, however, since the kinetic energy does 
not exceed the magnetic energy. As demonstrated by 
Galsgaard & Nordlund (1996), a magnetically dom- 
inated plasma is able to quickly dissipate structural 
complexity, independent of the value of the resistivity. 



Thus, a sub-Alfvenic field could not remain strongly 
tangled, and hence could not avoid detection. Using 
the same argument, we expect those cores where a 
strong (sub-Alfvenic) field has indeed been observed 
to have a relatively simple magnetic field structure. 

Model Ad2 can be used for the description of a 
typical molecular cloud with linear size of about 15 pc 
and (n) = 130 cm" 3 . At the time t = 1.0td yn , the 
relation B—n is then: 



B w 4.5/iG( : 
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(13) 
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Since the exponent is > 0.5, most of the dense cores 
are found to have magnetic pressure larger than ther- 
mal pressure, at early times. The lowering with time 
of the exponent in the B-n power law means that 
later on, the dominance of the magnetic field in dense 
cores tends to be reduced. Although the flow is 
random and approximately isotropic, the kinematic 
alignment of B with v makes dense cores accrete mass 
along B at an increased rate. Therefore, the accretion 
of mass around dense cores, embedded in a random 
flow with /3 rs 1, is such that, while magnetic pressure 
becomes dominant over thermal pressure during the 
initial phase of turbulent fragmentation, later on, in 
some cores, magnetic pressure can decreases to the 
level of thermal pressure, even in the absence of grav- 
ity, and on a time scale that is competitive with am- 
bipolar diffusion. This mechanism could be relevant 
for the process of star formation. 

A central point that we want to stress in this pa- 
per is that the fact of finding some cloud cores, with 
sub-Alfvenic velocity dispersions and with magnetic 
pressure larger than thermal pressure, does not nec- 
essarily mean that the dynamics of molecular clouds 
is dominated by MHD waves; those cores may be 
formed, in a few million years, in a super-sonic and 
super-Alfvenic flow, only marginally affected by the 
magnetic field (model Ad2). The dynamics becomes 
strongly affected by the magnetic field only in some 
very dense regions, on small scales, and preferentially 
during the first dynamical time. 

We stress that a key theoretical ingredient to the 
interpretation of OH Zeeman splitting data is the B- 
n relation. The fact that the exponent of the relation 
is > 0.5 for more than one dynamical time is the 
reason why sub-Alfvenic cores can be found in the 
experiment. 

An interesting question is whether the present re- 
sults, that indicate that kinetic energy dominates 
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over magnetic energy on the scale of cold clouds, 
could still be consistent with the well established 
rough equipartition of kinetic and magnetic energies 
on larger scales, e.g. comparable with the thickness 
of the galactic disc ( Boulares fc Cox, 1990| ). Here, 
we can only offer some preliminary remarks; a more 
quantitative answer requires larger scale and higher 
resolution numerical simulations, that go beyond our 
current numerical capabilities. 

However, we have made a rather interesting obser- 
vation already in the current simulations. It turns out 
that, while the power spectra of velocity and mag- 
netic field strength are nearly parallel, the spectra 
of kinetic and magnetic energy are not. Because of 
the intermittency of mass density p that develops in 
these supersonic flows, the power spectrum of p^v 
need not be parallel to the power spectrum of v. In 
our super-Alfvenic experiments there is a difference of 
slope between the power spectra of magnetic and ki- 
netic energy, amounting to about 0.5, with the kinetic 
energy spectrum being more shallow. 

Such a difference in slope opens up the interesting 
possibility that, over the 2-3 orders of magnitude that 
separate cold clouds from large scale disc structures, 
a difference of about an order of magnitude devel- 
ops between the kinetic and magnetic energy densi- 
ties. This is roughly consistent with the differences 
between kinetic and magnetic energy adopted in our 
super-Alfvenic models. More firm conclusion on this 
issue must await models with higher Mach number 
and higher numerical resolution. 

5. Conclusions 

In this work we have shown that: 

• Both super-Alfvenic and equipartition random 
super-sonic flows, are characterized by a de- 
cay time-scale approximately equal to one dy- 
namical time. One of the original motivations 
for the theoretical models of magnetized clouds 
with rough equipartition between kinetic, grav- 
itational, and magnetic energies was the belief 
that in such models the decay time-scale could 
be significantly longer than one dynamical time. 
Since this is now known to be incorrect (see 
also MacLow et al. 1998 and Stone, Ostriker & 
Gammie 1998), equipartition models of molec- 
ular clouds have lost their original motivation. 

• Random super-sonic and super-Alfvenic mo- 



tions produce a very intermittent probability 
distribution of magnetic field strength, with an 
exponential tail. Cores with field strength sev- 
eral times in excess of the mean field can be 
formed. 

• Random super-sonic and super-Alfvenic mo- 
tions also produce a very intermittent proba- 
bility distribution of gas density, in agreement 
with stellar extinction measurements, while 
equipartition models fail to reproduce the ob- 
servational result. 

• Synthetic spectra calculated from both the super- 
Alfvenic and the equipartition models are used 
to study the correlation between line width and 
gas column density. The equipartition models 
fail again to reproduce qualitatively the growth 
of line width with column density, while the 
super-Alfvenic model shows a correlation simi- 
lar to the observed one. 

• Dense cores, with magnetic pressure larger than 
thermal pressure, and velocity dispersions smaller 
than va , are found as the result of the evolution 
of supersonic and super-Alfvenic flows. 

• A power law statistical B-n relation is gen- 
erated by super-sonic and super-Alfvenic mo- 
tions, but not by MHD waves in the equipar- 
tition model. The exponent of the relation is 
> 0.5 for about one dynamical time, which al- 
lows for the existence of cores with va > Cs- 

• The exponent in the B-n relation of the super- 
Alfvenic model decreases with time, because 
magnetic field lines are stretched and partially 
aligned with the flow. The statistical impor- 
tance of the magnetic pressure in dense cores 
thus decreases with time, even in the absence 
of gravity and ambipolar diffusion. After about 
one dynamical time, the exponent of the upper 
envelope of the B-n relation, is approximately 
equal to 0.4, as in the observations. 

• The scatter in the B-n relation of the super- 
Alfvenic model is also consistent with the Zee- 
man splitting measurements, while the scatter 
in the equipartition model is inconsistent with 
most Zeeman splitting upper limits on the field 
strength, and even with some field detections. 
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In summary, even though the evidence presented in 
the present paper is far from complete, it appears that 
super- Alfvenic motions produces in a natural way the 
same type of statistical properties that characterize 
cold molecular clouds, while motions with equiparti- 
tion of kinetic and magnetic energies have properties 
that are harder to reconcile with the observations. 
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Figure and Table captions: 



Table 1: Initial parameters of the numerical experiments. 
Ma'- alfvenic rms Mach number; 
&o '■ magnetic field strength parameter; 
M: rms Mach number; 
ambi: ambipolar diffusion parameter; 



Table 2: Initial physical parameters of the models, obtained by scaling the numerical parameters of the exper- 
iments (Table 1) with the Larson's relations (Eqs. 7 and (8)). 

Fig. [j]: Time evolution of the ratio of the energy decaying time divided by the instantaneous dynamical time. 
Separate plots are for the total energy E to t, the magnetic energy E m , and the kinetic energy E^. 

Fig. [2|: Same as in Fig. [I], but for models with ambipolar diffusion. 

Fig. ^: Distribution of magnetic field strength in the super-Alfvenic model Ad2 (thick line), and in the equipar- 
tition model Bdl (thin lines). The dashed vertical lines mark the mean field strength in the two models. 

Fig. ^: Histograms of extinction (left panels) and plots of dispersion of extinction in cells versus the mean cell 
extinction (right panels). The top panels show the result from the super-Alfvenic model Ad2, while the middle 
panels are from the equipartition model Bdl. The bottom panels are Lada ct al. (1994)| observational data for the 
cloud IC5146. 



Fig. |5f Histograms of integrated antenna temperature of synthetic J=1^0 13 CO spectra, from the super- 
Alfvenic model Ad2 (thick line) and for the equipartition model Bdl (thin line). Model Ad2 shows a very inter- 
mittent approximately exponential distribution. 

Fig. ^: J=1^0 13 CO equivalent width versus J=1^0 13 CO integrated antenna temperature, from the super- 
Alfvenic model Ad2 (upper panel) and for the equipartition model Bdl (lower panel). 

Fig. [7]: The B — n relation: observations and theoretical models. The thick contour lines are from the the 
super-Alfvenic model Ad2, and the thin contour lines from the equipartition model Bdl. 

Fig. |^: Three dimensional visualizations of the density field in the super-Alfvenic model Ad2, obtained by the 
superposition of voxel projection and density isosurfaces. The left panel corresponds to t = 0.5td yn ,o, and the right 
panel t = lM dvn ,o- 

Fig. ^|: Same as in Fig. ^|but for the equipartition model Bdl. The left panel corresponds to t = l.3td yn ,o, an d 
the right panel t = 2.8td yn ,o- 

Fig. [u]: Histo grams of the cosine of the angle between v and B. In experiment Al (left panel) there is a partial 
alignment, while in experiment Bl (right panel) the two fields are mainly perpendicular to each other. 
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Model name 


<7„/(km/s) 


Z/pc 


tW(ioV) 


<n}/(100 H 2 cm- 3 ) 
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